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Abstract. For inviscid fluid flow in any n-dimensional Riemannian manifold, new con- 
served vorticity integrals generalizing helicity, enstrophy, and entropy circulation are derived 
for lower-dimensional surfaces that move along fluid streamlines. Conditions are determined 
for which the integrals yield constants of motion for the fluid. In the case when an inviscid 
. fluid is isentropic, these new constants of motion generalize Kelvin's circulation theorem 

from closed loops to closed surfaces of any dimension. 
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1. Introduction 

Vorticity conservation laws have long been of interest in the study of inviscid fluid flow 
(e.g. [HOE]). Their mathematical formulation in multi-dimensions is given by an integral 
continuity equation 

(/ f TdV = - [ X-dA (1.1) 

'V(t) JdV(t) 



dt 



for a conserved density T and a spatial flux X that have an essential dependence on the curl 
of the fluid velocity, as defined on domains V(t) that are transported along the streamlines 
of the fluid governed by Euler's equations. Physically speaking, such conservation laws 
express basic rotational properties of the fluid in a reference frame that moves with the 
fluid flow jl]. To-date the only known vorticity conservation laws consist of the general 
enstrophy integral (5) H |7] which exists in two (and higher-even) spatial dimensions, the 
helicity integral jH El [7] which exists in three (and higher-odd) spatial dimensions, and the 
entropy circulation integrals [91 [10] which exist in two as well as three (and higher) spatial 
dimensions. A general formulation of these conservation laws, applicable to incompressible 
as well as compressible inviscid fluids, is hard to find in the literature. 

The Eulerian fluid equations in M n are given in terms of the velocity u, the mass density 
p, the entropy density S, and the pressure p by 

u t + u ■ Vu = -p~ l Vp, (1.2) 
p t + V ■ (pu) = 0, (1.3) 
S t + u-VS = 0. (1.4) 
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In the case of compressible fluid flow, p is specified by an equation of state in terms of a 
function p = P(p, S), whereas in the case of incompressible fluid flow, p and S are constant 
while p satisfies the Laplacian equation —p^Ap = (Vtt) ■ (Vu). 

For inviscid fluid flow in R 2 , enstrophy and entropy circulation are defined by the respective 
two-dimensional integrals 

/ p~ l w 2 d 2 x (1.5) 

JV(t) 

and 



Szu d x = <b Su-ds + / u-*VSdx (1.6) 
v(t) Jav(t) Jv(t) 

where w = V ■ *u is the vorticity scalar, given in terms of the infinitesimal rotation operator 
*. Enstrophy physically measures the bulk rotation of the fluid, while entropy circulation 
measures the flux of entropy around closed loops when the streamlines are aligned with the 
entropy gradient in the fluid. The more general enstrophy integral 

f{w/p)pd 2 x = I j{wjp)wd 2 x (1.7) 
V(t) JV(t) 

is connected with circulation properties of the fluid, when the function / is odd (/ is even), 
and with bulk rotation moments of the fluid, when the function / is even (/ is odd). 

In contrast, for inviscid fluid flow in M 3 , helicity and entropy circulation are respectively 
defined by the three-dimensional integrals 

u-Sd 3 x (1.8) 

V(t) 



and 



VS-$d 3 x= (b S$-dA= <b (ux VS) ■ dA (1.9) 

V(t) JdV(t) JdV(t) 



involving the vorticity vector •& = V x u, given in terms of the cross-product operator. 
Physically, helicity measures the net rotation of the fluid around the directions of streamlines, 
while entropy circulation detects the net alignment between the streamlines and the entropy 
gradient. A more general integral 

I f(VS ■$/p)pd 3 x= [ f(VS -^/p)VS -{}d 3 x (1.10) 

JV(t) JV(t) 

measures entropy circulation properties of the fluid, when the function / is odd (/ is even), 
and bulk alignment between the streamlines and the entropy gradient in the fluid, when the 
function / is even (/ is odd). 

The helicity integral (jl.8p and the general enstrophy integral (ll.7j) are conserved whenever 
the fluid flow is isentropic (i.e. S is constant throughout the fluid domain) and either 
incompressible (i.e. p is constant throughout the fluid domain) or compressible with a 
barotropic equation of state in which the pressure p is a function only of the density p, 
whereas the entropy circulation integrals (II. 6ft and (ll.lOj) are non-trivially conserved only 
when the fluid flow is non-isentropic (i.e. S is constant only along streamlines) as well as 
compressible (so that p is non-constant in the fluid domain). All of these conserved vorticity 
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integrals have a mathematical origin as Casimir invariants of the Hamiltonian formulation 
PU CLl El Q31 E] for the underlying fluid equations (0]) - (Oj) . 

In the present paper, new conserved vorticity integrals generalizing helicity, enstrophy, 
and entropy circulation are derived for p-dimensional surfaces that move along streamlines 
of inviscid fluid flow in any n-dimensional Riemannian manifold with n > p > 1. In the 
case of odd-dimensional surfaces, the helicity integral yields a constant of motion if the 
surface is boundaryless or satisfies a zero-vorticity boundary condition, while the entropy 
circulation integral yields a constant of motion whether the surface is boundaryless or has a 
boundary. Similarly, in the case of even- dimensional surfaces, the entropy circulation integral 
yields a constant of motion if the surface is boundaryless or satisfies a vorticity boundary 
condition involving the entropy gradient vector, while the enstrophy integral yields a constant 
of motion whether the surface is boundaryless or has a boundary. As a by-product of these 
results, a simple proof of Kelvin's circulation theorem is obtained for isentropic fluid flow 
in all dimensions n > 1, and its relationship to the new helicity and enstrophy integrals for 
two-dimensional moving surfaces is explained. 



2. Inviscid fluid equations in Riemannian manifolds 

Consider an n-dimensional manifold M with a Riemannian metric g. Let V be the metric- 
compatible covariant derivative determined by Vg = 0, and let e g be the metric- normalized 
volume tensor determined by Ve g = and g(e g ,e g ) = n\. Write 9 V and div for the vector 
derivative operator and the covariant divergence operator defined by £JV = g(£, 9 V) and 
g( 9 V, £) = div£ holding for an arbitrary vector field £ on M. These operators are the natural 
covariant counterparts of the gradient V and divergence V- operators in R n . 

In this tensorial notation, the covariant generalization of Euler's equation (11. 2p from M. n 
to M is given by 

u t + {u\V)u = -p- l9 Vp (2.1) 

where u is the fluid velocity vector on M. Similarly the covariant equations for the fluid 
pressure p, mass density p, and entropy density SonM are given by 

p = P(p,S), (2.2) 
p t + div(pu) = 0, (2.3) 
S t + u\VS = 0, (2.4) 

in the case of compressible fluid flow, or by 

p = const., S = const., (2.5) 

- p- 1 A g p = g( 9 \7u, 9 Vu) + R(u,u) (2.6) 

in the case of incompressible fluid flow, where A g = g( 9 V, 9 V) = 5 VJ V denotes the scalar 
Laplacian, and R denotes the Ricci tensor of g. In both cases, the curl of the velocity is an 
antisymmetric tensor on M 

u = 9 V Au (2.7) 

obeying the transport equation 



u t + 9 VAg(u,u) = p- 29 VpA 9 Vp, } VAw = 0. 
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(2.8) 



It will be very advantageous to rewrite the fluid equations (I2.1l) - fl2.8p by means of differ- 
ential forms combined with the material (advective) derivative 

T> t = D t + C u . (2.9) 

Here D t denotes dt acting as a total time derivative, and C u denotes the Lie derivative acting 
on p-forms a on M by C u oc = u\ da + d(wj ct). Note that the exterior derivative d commutes 
with the Lie derivative C u and acts as a total spatial derivative. Now let u be the velocity 
1-form and let us be the curl 2-form defined by the respective duals of u and u with respect 
to the metric g on M, namely £ju = g(£,,u) and (\u: = g((,u:) for an arbitrary vector field 
£ and an arbitrary antisymmetric tensor field ( on M. The Lie derivatives of u and u) along 
streamlines are, respectively, 

C u u = u\uj + d(u\ u), C u u = d(u\uj) 

where 

u\u> = (u\ V)u — d(iiiju). 
Then the equations for S, p, u, and uj have the elegant transport formulation 

® t S = 0, 
£) t p = -pdiv g u, 

S t u = d(||u|g — e — p^p) + esdS, du = u>, 
D t uj = de s A dS, du: = 0, 

with diVgVL = 9 Vju = divu and |u|^ = u\u = g(u,u), where 

e = J p- 2 pdp (2.16) 

is the thermodynamic energy of the fluid. In particular, note e = J p~ 2 P(p, S)dp is a 
specified function of p and S for compressible flow, whereas e = (since dp = 0) vanishes 
for incompressible flow. 



(2.10) 
(2.11) 

(2.12) 
(2.13) 
(2.14) 
(2.15) 



3. Moving surfaces and conserved integrals 

Now consider any smooth orientable p-dimensional submanifold S(t) G M that is trans- 
ported along streamlines in the fluid. In particular, in local coordinates, each point x l G S(t) 
obeys dx % /dt = C u x l = u\Vx % (i = 1, . . . , n). A conserved integral on the moving surface 
S(t) consists of an integral continuity equation given by 




for a p-form density o: and ap - 1-form flux /3 that are functions of the space and time 
coordinates x l , t, and the fluid variables S, p, u, uj (and possibly their spatial derivatives) 
subject to the fluid equations (I2.12p - fl2.14p . 

Integrals of this form (13. ip will define a constant of motion on the moving surface S{t) for 
all formal solutions of the fluid equations provided that the flux integral on dS(t) vanishes 
identically. If S(t) is boundaryless then every conserved integral (13.1 p yields a constant 
of motion. Alternatively, if S(t) has a boundary then a conserved integral (13.1 ft yields a 
constant of motion only when f3 = (modulo an exact p—1 form). Note that S(t) will be a 
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moving domain V(t) if p = n, in which case the conserved integral (13. ip coincides with the 
continuity equation 

— ! TdV = - [ g(X, v)dA (3.2) 
dt Jv(t) Jav{t) 

given by 

a = eT, (3 = -X\e (3.3) 

where dV = e g is the volume n-form dual to the tensor e g , and dA = z>Je 3 is the hypersurface 
area n — 1-form in terms of the unit-normal vector v. In this situation the conservation law 
(13 .4p reduces to a space-time divergence D t T + DivX = where Div denotes div acting as 
a total derivative. 

Proposition 1. A p-form density ct and ap — 1-form flux f3 yield a conserved integral ( 13. 1ft 
on all p- dimensional moving surfaces S(t) C M iff 

D t a = d/3 (3.4) 

holds for all formal solutions of the fluid equations. 

Proof: Let <\> t be the diffeomorphism of M defined by the streamlines of the vector field 
u\ with 0o = id. The moving surface S(t) can then be viewed as the flow of a fixed 
surface S(0) in M. Under this flow, the pullback of a \ s ^ is given by 0*a| 5 ^, where 
d 

Jt [ 



f t ct) = D t cx. + £ u ot = ®tOL. Hence we have 



d 



a 

>S{t) 

and thus the integral continuity equation ( 13. ip becomes 



t=o JS{o) dt js(o) 



® t cx (3.5) 



0=/ V t a- (3= (J) t a-dp). (3.6) 

JS(0) JdS(0) Js(0) 

The vanishing of this integral for an arbitrary p-dimensional surface S(0) is equivalent to 
the conservation law (13.41) . ■ 

Conserved vorticity integrals in which the density a and the flux (3 have an essential 
dependence on the curl 2-form a; will now be derived from the transport equations (I2.12p - 
( I2.15p . first for isentropic fluid flow and then for non-isentropic fluid flow. 

4. Helicity on moving surfaces 

Three-dimensional helicity (I1.8P can be naturally generalized to all odd dimensions n > 3 
in terms of the vorticity vector [7J [3] 

= = t g \u} {n - 1)/2 (4.1) 

(n— 1)/2 times 

defined on the n-dimensional Riemannian manifold M. Here * is the Hodge dual operator, 
acting by contraction with respect to the volume tensor e g , and powers of u) denote exterior 
products. For a moving domain V(t) C M, the n-dimensional helicity integral for isentropic 
fluid flow in M is given by a continuity equation (13. 2p holding for the conserved density T = 
0Ju and the spatial flux X = — cr$ with o = ||u|^ — e — p/p, where e is the thermodynamic 
energy of the fluid. A simpler formulation is obtained by using the corresponding differential 
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forms (I3.3p : a = u A uj^ 1 )/ 2 and j3 = cru)^ 1 '' 2 . The resulting conserved helicity integral 
( 13. ip has a natural extension from moving domains in M to odd-dimensional moving surfaces 
in M, whether n is odd or even. 

Theorem 2. For isentropic inviscid fluid flow in any Riemannian manifold M with dimen- 
sion n > 1, the 2q + 1-form density a = u Aw 5 and the 2q-form flux (3 = (^\u\ 2 + e — p/ p)u) q 
satisfy the continuity equation 

^ / uAw'= / (i|ug + e - p/p)w«, < Q < (n - l)/2 (4.2) 

at JS(t) JdS(t) 

yielding a conserved helicity integral on any smooth orientable submanifold S(t), with odd 
dimension 2q + 1, transported along streamlines in the fluid. 

Proof of conservation: Isentropic fluid flow has dS = in M, so thus the velocity and curl 
equations (I2.14l) - ()2.15p reduce to 

D t u: = 0, D 4 u = dcr, a = ||u|J - e - p~*p (4.3) 

with 

duj = 0, du = u). (4.4) 

Since ®tu is an exact 1-form and a; is a closed 2-form annihilated by D t , we have T>t(uAu} q ) = 
da A ua q = d(cruj q ), which establishes (jOl) . ■ 

When the surface S is a closed curve (i.e. dS=0 with q — 0), the helicity integral (14.21) 
reduces to the circulation integral 

*L u=0 (45) 

providing a generalization of Kelvin's circulation theorem to isentropic fluid flow (compress- 
ible or incompressible) in Riemannian manifolds of any dimension n > 1. 



5. Enstrophy on moving surfaces 

Two-dimensional enstrophy ( jl.7p has a natural extension to all even dimensions n > 2 in 
terms of the vorticity scalar [3 [3] 

zj = *( pA- - -Au>) = e q \uj n/2 (5.1) 

n/2 times 

defined on the n-dimensional Riemannian manifold M in terms of the Hodge dual operator 
*. Here powers of u> again denote exterior products. For a moving domain V(t) C M in 
isentropic fluid flow, the n-dimensional enstrophy integral is given by a continuity equation 
( 13. 2 p holding for the conserved density T = f(w/p)w with no spatial flux X = 0, where / 
is an arbitrary non-constant function. 

An equivalent formulation uses the corresponding differential forms ( 13. 3p : a. = 
f(w/p)uj n l 2 and (3 — 0. The resulting conserved enstrophy integral (13. ip can be gener- 
alized from moving domains in M to even-dimensional moving surfaces in M. 
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Theorem 3. For isentropic inviscid fluid flow in any Riemannian manifold M with even 
dimension n > 2, the 2q-form density ct = fuj q depending on an arbitrary non-constant 
function fiwjp) satisfies the continuity equation 



d 
dt 



fu q = 0, 1 < q < n/2 (5.2) 

s{t) 



yielding a conserved enstrophy integral on any smooth orientable submanifold S (t) , with even 
dimension 2q, transported along streamlines in the fluid. 

Proof of conservation: The curl equation (14. 3[) holding for isentropic fluid flow implies that 
®t{fv q ) = ®t(^/p)f'u q , where w = e g \uj n l 2 . Then using 

D t w = (D t e g )\uj n/2 , D t e g = £ u e g = -(divu)e g , D t p = -{divu)p, (5.3) 

we have 

©tOT 1 *,) = (5.4) 

whence D t (zu/p) = 0. This implies D t (fuj q ) = 0, establishing AO]) . ■ 
The enstrophy integral (15.21) can be alternatively expressed as 

±[ fu q = -^-f f&{w/p) A u A w*" 1 + ± I /uA^ = (5.5) 

at JS(t) at JS(t) at JdS(t) 



through the properties (14.41) of u), whether the surface S has a boundary or is boundaryless. 
Consequently, if the function / is constant, then the enstrophy will reduce to the helicity on 
the boundary of S or will otherwise vanish when S has no boundary. 



6. Entropy circulation on moving surfaces 

Similarly to enstrophy and helicity, two-dimensional entropy circulation (11.61) and three- 
dimensional entropy circulation (11.101) have a straightforward extension respectively to all 
even dimensions n > 2 and all odd dimensions n > 3. The entropy circulation integral for 
a moving domain V(t) C M in non-isentropic fluid flow is given by a continuity equation 
(13. 2p that holds in the even- dimensional case for the conserved density T = f(S)w and the 
spatial flux X = —f(S)es(dw/du:)\dS, where e is the thermodynamic energy of the fluid, 
and that holds in the odd- dimensional case for the conserved density T = f(r/p)r with 
vanishing spatial flux X = 0, where 

r = *(ds Aw (n ~ 1)/2 ) = $\ds (6.i) 

defines the entropy circulation scalar on the Riemannian manifold M in terms of the vorticity 
vector *&. The corresponding differential forms (I3.3P are given by a = f{S)u n / 2 and (3 = 
^nf(S)esdS Aa/ n ~ 2 )/ 2 when the dimension n is even, and by ex. = f(r/p)iv^ n ~ 1 ^ 2 and (3 — 
when the dimension n is odd. In each case, there is a natural generalization of the resulting 
conserved entropy circulation integral (13. ip to moving surfaces in M. 
The generalization is presented first for even- dimensional surfaces. 

Theorem 4. For non-isentropic inviscid fluid flow in any Riemannian manifold M with 
even or odd dimension n > 2, the 2q-form density a = fuj q and the 2q — 1-form flux (3 = 
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qesfdS A uj q 1 depending on an arbitrary non- constant function f(S) satisfy the continuity 
equation 

-I / fu,* = q<[ esfdSAu*- 1 , l<q<n/2 (6.2) 

at JS(t) JdS{t) 

yielding a conserved entropy circulation integral on any smooth orientable submanifold S(t) , 
with even dimension 2q, transported along streamlines in the fluid. 

Proof of conservation: From the entropy equation ( 12.12ft and the curl equation ( 12.15ft 
holding for non-isentropic fluid flow, we have Q t (f<jJ q ) = qfdes A dS A uj^ 1 = qd(esfdS A 
oj^ 1 ) since the 2-form uj is closed. This establishes (I6.2p . ■ 

When the surface S is boundaryless, the entropy circulation integral (I6.2p becomes 

!Li fu q = -^-l f'dS A\iAut 9 - 1 = (6.3) 
dt Js(t) dt Js{t) 

from the properties (14.4ft of u>. As a consequence, if the gradient of the entropy S is aligned 
with the fluid streamlines on the surface S (which will occur whenever the fluid flow is 
isentropic in S C M), then the entropy circulation (16. 3p will vanish due to dS A u = 0. 
The analogous generalization for odd- dimensional surfaces is presented next. 

Theorem 5. For non-isentropic inviscid fluid flow in any Riemannian manifold M with 
odd dimension n > 3, the 2q + 1-form density a = fdS A uj q depending on an arbitrary 
non- constant function f(r/p) satisfies the continuity equation 

[ fdSAu q = 0, l<?<(n-l)/2 (6.4) 
Js(t) 



dt 



yielding a conserved entropy circulation integral on any smooth orientable submanifold S(t) , 
with odd dimension 2q + 1, transported along streamlines in the fluid. This integral remains 
conserved if f is generalized to be a function of both Fjp and S. 

Proof of conservation: The curl equation ( 12.15ft combined with the entropy equation ( 12.12ft 
yield T> t {dSAoj q ) = and ® t f = f t D t r = f t ® t {p~ l e g )\ (dSAu/™- 1 )/ 2 ) with f = Fjp where 
F = e g \ (dS A u^- 1 ^ 2 ). Then equation (JED implies D t f = and hence D t (fdS A u q ) = 0, 
which establishes (16. 4p . ■ 

In contrast to the even- dimensional integral (16. 2p . here the odd- dimensional entropy cir- 
culation integral (16 .4p can be expressed as 

^ / fdS A = ~ [ f'Sd{r/p) A + ^ / fS<J = (6.5) 

at JS(t) at JS(t) at JdS{t) 

whether the surface S has a boundary or is boundaryless. If the function / is constant, 
then this integral (I6.5P will reduce to the entropy circulation on the boundary of S or will 
otherwise vanish when S has no boundary. 



7. New constants of motion 

As a preliminary step, for any s- dimensional smooth orientable surface S C M, let 
{ e ([i)}fj,=i,...,n-s be an arbitrary orthonormal basis for the normal space of S in T X M at each 
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point of S. Then the volume forms of the surface S and its boundary dS can be defined by 
the projections 

dV s = e(S) = (e {1) A • • ■ A e (n _ s) )J e g (7. 1) 

and 

dA s = e(5) = e(dS) = (z> 5 A e ( i) A ■ ■ ■ A e (n _ s) )J e 9 (7.2) 

where i>s denotes the unit normal of dS in S and e g is the volume n-form of M. Additionally, 
the dual volume tensors e(S) for S and e(dS) for 95 can be defined in terms of the volume 
tensor e g for M by 

e{S) A e ( ij A • ■ • A e (n _ s) = e p , e(0S) A z> 5 = e(«S). (7.3) 

Now, if S(t) is an odd- dimensional surface transported along fluid streamlines in M, the 
helicity integral f)4.2p for isentropic fluid flow and the entropy circulation integral (16. 4p for 
non-isentropic fluid flow can be formulated respectively as 

' [ # s \u s dVs=<( {\\u\ 2 g -e-p/p)g{'ds,i>s)dAs (7.4) 

JS(t) JdS(t) 



dt 



and 



Here 



^/ f(S,r/p)# s \dsSdV s = 0. (7.5) 
at Js(t) 



$ s = 6(S)\(u\ s )i (7.6) 

defines the vorticity vector of the surface S with q = |(dim(«S) — 1), while 11,5 and d s S 
denote the respective projections of the velocity 1-form u and the gradient 1-form d5* into 
the cotangent space of S, and r is the entropy circulation scalar (16. lft of the manifold M 
when the dimension n is odd. 

Similarly, if S(t) is an even-dimensional surface transported along fluid streamlines in M, 
the enstrophy integral (I5.2p for isentropic fluid flow and the entropy circulation integral (16.21) 
for non-isentropic fluid flow have the respective formulations 

^ / f{w/p)w s dV s = (7.7) 
at Jstt) 



and 



where 



^ / f{S)w s dV s = q(f e s f(S)g($ ds , V 9S S) dA s (7 1 

at Js(t) JdS(t) 



w s = e(S)\(uj\ s y (7.9) 

defines the vorticity scalar of the surface S with q = |dim(«S), while Vq$ denotes the 
projection of the gradient 9 V into the tangent space of dS, and w is the vorticity scalar 
(15. ip of the manifold M when the dimension n is even. 

These formulations (I7.4p — (l7.9p are useful for investigating the constants of motion that 
arise under various conditions for moving surfaces S(t) in isentropic and non-isentropic fluid 
flow in M as follows. 
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7.1. Closed moving surfaces spanned by a moving hypersurface. Consider any even- 
dimensional surface £(£) spanning a closed orientable surface S(t) transported along the fluid 
streamlines in M, with diniiS = dimE — 1, <9£ = S and dS = 0. Then Stokes' theorem can 
be applied to convert the conserved helicity integral (17. 4p into the equivalent form 

-I #s\uisdVs = ^-f zu s dVx = (7.10) 

dt JS{t) dt JT.(t) 

yielding a scalar vorticity constant of motion on S(t) for n-dimensional isentropic (compress- 
ible or incompressible) fluid flow. This new constant of motion (I7.10P provides a higher- 
dimensional version of Kelvin's circulation (14. 5p . which has the analogous formulation as a 
vorticity constant of motion 

AC AC 

' Us ds = — zux dV^ = (7.11) 



dt J c dt j S (t) 

on 2-surfaces S(t) that span a closed curve C(t) transported along the fluid streamlines in 
any even or odd dimension n > 2. (Note here us = e(C)Ju denotes the scalar projection of 
u along the curve C, and ds = e(C) is the arclength 1-form of C.) When the dimension n is 
even, these constants of motion (17.101) and (17. lip are special cases of the enstrophy constant 
of motion (17. 7p on E(t). Therefore, the helicity and enstrophy integrals can be unified in the 
form of a generalized scalar vorticity constant of motion 

^ / f{w/p)wv dVx = (7.12) 

JY.(t) 

in which / is a function of w/ p when n is even or / is a constant when n is odd. 

In contrast, for n-dimensional non-isentropic (compressible) fluid flow, the scalar vorticity 
integral (17. 12[) is no longer conserved but instead is replaced by the conserved entropy circu- 
lation integrals (I7.5P and (17. 8p . Given any odd-dimensional surface S(t) spanning a closed 
orientable surface S(t) transported along the fluid streamlines in M, with dim S = dim S — 1, 
<9£ = S and dS = 0, the integral (17.81) can be converted through Stokes' theorem into a 
special case of the integral (17. 5p on E(t) when the dimension n is odd, 

f{S)w s dV s = ^ + f f'(S)$x\dSdVx = 0. (7.13) 

As a result, both entropy circulation integrals (17. 5p and (17.81) have a unified formulation 

^/ f(S,r/p)^\dSdV s = (7.14) 

at JT,(t) 

yielding a new vorticity constant of motion on S(t) for n-dimensional non-isentropic (com- 
pressible) fluid flow. This constant of motion (I7.14p measures the alignment between the 
entropy gradient and the vorticity vector of the fluid on £(i), weighted by a function / of S 
when n is even or a function / of both S and F j p when n is odd. 

7.2. Non-closed moving surfaces and vorticity boundary conditions. Consider an 
orientable non-closed, but otherwise arbitrary, surface S(t) transported along the fluid 
streamlines in M. Then the entropy circulation integral (17.51) when S is odd- dimensional 
and the enstrophy integral (17. 7p when S is even-dimensional are constants of motion respec- 
tively for non-isentropic fluid flow and isentropic fluid flow in M. In contrast, the helicity 
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integral (17. 4p when S is odd- dimensional and the entropy circulation integral (17. 8p when S 
is even-dimensional are constants of motion if and only if the flux integrals on the boundary 
dS(t) vanish respectively for isentropic fluid flow and non-isentropic fluid flow in M. 

The helicity flux integral (I7.4p clearly vanishes if g($s, ^s)\ds — holds on the boundary 
surface dS(t). This condition states that the vorticity vector $5 at dS(t) must lie in the 
tangent space of dS(t). A useful equivalent formulation is that the vorticity scalar of the 
boundary surface dS(t) must vanish everywhere on dS(t), 

Was = (7.15) 

since g($s, As) las — ro as holds due to the identity (17. 3p relating the volume tensors of S{t) 
and dS(t). 

Likewise, the entropy circulation flux integral (17. 8p vanishes whenever g(^as,^asS)\ds = 
0, so that the vorticity vector figs must be orthogonal to the entropy gradient at the boundary 
surface dS(t). This condition has the equivalent formulation 

$ ds \d ds S = (7.16) 

where dg$ denotes the projection of the exterior derivative into the cotangent space of dS. 

As a consequence of the following result, these two flux conditions (I7.15P and (I7.16P are 
readily seen to be preserved under transport along fluid streamlines. 

Lemma 6. For isentropic fluid flow in M , 

®t™as = -UJ ds dtv 9s u, (7.17) 

and for non-isentropic fluid flow in M , 

®t(ti8s\d 9S S) = -(# ds \d dS S)dw dS u, (7.18) 

where divgs denotes the divergence operator projected into the boundary surface dS. 

Proof: View S(t) as the flow of a fixed surface S(0) under the diffeomorphism <p t of 
M generated by the streamlines u\ t with O = id- The pullback of the boundary vol- 
ume form dAs = e(dS(t)) of dS(t) yields — (fitdAgl = C u dA s \ = (div 9s u)dA s on 

the boundary surface dS(0). This implies C u e(dS) = {divQsu)e(dS) on dS, and hence 
the dual volume tensor e(dS) obeys C u e(dS) = —(divQsu)e(dS). Now let q = dim dS. 
From teas = e(dS)\uj q ^ 2 , we have Dt^as = C u e(dS)\uj q ^ 2 = —{divQ S u)e(dS)\u; q / 2 since 
T) t ijJ = holds for isentropic (compressible or incompressible) fluid flow. Similarly, from 
#as\dasS = e(dS)\ (dS A o;^ 1 )/ 2 ), we get T>0 9s \d 9s S) = C u e(dS)\ (dS A a;^ 1 )/ 2 ) = 
— (div 9 su)e(dS)\ (dS A uA -1 )/ 2 ) since T) t (oj A dS) = holds for non-isentropic fluid flow. ■ 
Therefore, when an orientable non-closed odd-dimensional surface S(t) satisfies the vortic- 
ity boundary condition (I7.15P transported along fluid streamlines in M, the helicity integral 
(I7.4p yields a constant of motion 



d 
dt 



$ s \n s dVs = (7.19) 

S(t) 



for isentropic fluid flow in M. Similarly, when an orientable non-closed even-dimensional sur- 
face S(t) satisfies the vorticity boundary condition (I7.16P transported along fluid streamlines 
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in M, the entropy circulation integral (17. 8p yields a constant of motion 



^ / f{S)w s dV s = 

at J S(t) 



(7.20) 



for non-isentropic fluid flow in M. 



8. Concluding remarks 



There are several interesting directions for future work. 

First, the physical meanings of helicity, enstrophy, circulation and entropy circulation are 
fairly well understood when M = R 3 ,R 2 . In higher dimensions, what is the precise physical 
content of the new constants of motion? This could be elucidated by evaluating the integrals 
(17.1 Op and ( I7.14p in M = R n for some physically relevant fluid configurations as well as for 
analytically interesting exact solutions of the fluid equations. 

Second, the helicity integral for a moving domain V(t) in M = R 3 is well-known to 
equal the average linking number of the integral curves of the vorticity vector i9 [5] . What 
is the relation between the new helicity/circulation constants of motion (IT. lOf) and ( 17. 19ft 
for moving surfaces in M = 1" with n > 3 and the topological linking of vorticity-lines? 
Likewise, is there an interpretation of the new entropy circulation constants of motion (I7.14p 
and (17.201) in terms of the topological structure of vorticity-lines and entropy gradient-lines 
for moving surfaces in M = R n ? 

Third, what information do these new constants of motion contain when the manifold M 
or the moving surface have nontrivial homology? 

Finally, another interesting open question is whether the fluid equations (I2.ip -( 12TS|) admit 
any additional conserved integrals that yield constants of motion on moving surfaces in n > 1 
dimensions, other than the conserved mass integral 



for n-dimensional domains V(t) transported along fluid streamlines in a Riemannian manifold 



It is worthwhile to write the new conserved integrals (14.21) . (15.21) . (16.21) and (16.41) in terms 
of the physical fluid variables u % and w' 3 ' = 2V' l, u^ given by local coordinates x % on M. Note 
the coordinate formulation of the transport equations (12.141) and (I2.15P for these variables 
is simply 




V(t) 



pdV = 



(8.1) 



M. 



Appendix A. Coordinate formulation 



u\ + u 3 Vju 1 = -p- 1 V> 

ul j + u k V k uj i3 - 2co k [l a 3]k = p- 2 V [l P V ]] p, V [k co ij] = 

where a l ° = 2V^ l u^ is the symmetric derivative of u 1 . 

Now consider isentropic fluid flow, i.e. S = const, in M, with 

p = P(p), p t + V l (pu l ) = 

when the fluid is compressible, or 



(Al) 
(A.2) 



(A.3) 



p = const., -p x V l Vjp = ViU j VjU l + RijU l u j 
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(A.4) 



when the fluid is incompressible, where Rij is the Ricci tensor and ViU^V jU l = \{pijG % i — 
UijU^) is the difference of the norms of the symmetric derivative <7 Jjf = 2W^u^ and the curl 
of u l . Then the helicity integral ( 14. 2ft for any orientable odd- dimensional surface S(t) 
transported along fluid streamlines in M is given by 

4 / QijU^s dV s = I {\9kiu k u l - e - p/ p)g ij d i s v j s dA s (A.5) 

at JS(t) JdS(t) 

in terms of the vorticity vector $ s of the surface S and the projection u l s of the velocity 
vector u % into the tangent space of S. Similarly, for any orientable even- dimensional surface 
S(t) transported along fluid streamlines in M, the enstrophy integral (15. 2p is given by 

4 / f(zu/p)zu s dV s = (A.6) 

at JS{t) 

in terms of the vorticity scalar rug of the surface S, where w is the vorticity scalar (15. ip of 
the manifold M when the dimension n is even. 

Next consider non-isentropic fluid flow, i.e. S ^ const, in M, with 

p = P(p,S), Pt + V l (pu l ) = 0, S t + uW i S = 0. (A.7) 

Then for any orientable odd- dimensional surface S(t) transported along fluid streamlines in 
M, the entropy circulation integral (16 ,4p is given by 

4 / f(r/p) 9ij ^V j s S dV s = 0, (A.8) 

at JS{i) 

where r is the entropy circulation scalar (16. ip of the manifold M when the dimension n is 
odd, while for any orientable even-dimensional surface S(t), the entropy circulation integral 
(16. 2p is given by 

4 / f(S)w s dV s = q [ e s f(S) gij 4 s V j dS S dA s (A.9) 

at JS(t) JdS(t) 

in terms of the projections of the gradient V* into the tangent spaces of S and dS. 
In all of these integrals, 

dV s = e h ... js (S)dx jl A • • • A dx js (A.10) 

is the volume element of the surface S, and 

dA s = n i e ijl ... js _ l (S)dx jl A • • • A dx^- 1 (A.ll) 

is the volume element of the boundary dS, with s = diniiS, where 6j v ..j s (S) denotes the 
metric-normalized volume form of S and h l denotes the unit normal of OS in S. 
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